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Abstract — Rational Bezier functions are widely used as 
mapping functions in surface reparameterization, finite element 
analysis, image warping and morphing. The injectivity (one-to- 
one property) of a mapping function is typically necessary for 
these applications. Toric Bezier patches are generalizations of 
classical patches (triangular, tensor product) which are defined 
on the convex hull of a set of integer lattice points. We give 
a geometric condition on the control points that we show is 
equivalent to the injectivity of every 2D toric Bezier patch with 
those control points for all possible choices of weights. This 
condition refines that of Craciun, et al., which only implied 
injectivity on the interior of a patch. 
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I. Introduction 

Mapping functions play an important role in computer 
graphics, computer aided geometric design (CAGD), finite 
element analysis (FEA) and some related areas. The injec- 
tivity of mapping functions, that is, the absence of self- 
intersection, is crucial in image warping and morphing ifTTI . 
free form deformation [1], surface reparameterization, and so 
on. Many authors have investigated conditions which imply 
injectivity. Goodman and Unsworth |7| proposed a sufficient 
condition for the injectivity of a 2D Bezier function. For 
the control points of a m x n tensor product patch, their 
condition involves 2m(m+l)+2n(n+l) linear inequalities. 
For image morphing, Choi and Lee [ 1 1 presented a sufficient 
condition for the injectivity of 2D and 3D uniform cubic B- 
spline functions. Their condition provides a single bound 
for the displacements of control points that guarantees the 
injectivity of the cubic B-spline function. Floater [6| studies 
a sufficient condition for injectivity of convex combination 
mappings over triangulations. 

Fig. Q] displays rational plane cubic Bezier curves with 
their control polygons (bold lines). The curve in Fig 
has no points of self-intersection. The curve in Fig. 
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has one point of self-intersection, which may be removed 
by varying the weights as shown in Fig. |l(c)| The control 
polygon of the first curve is in convex position, so there 
are no positive weights for which the resulting Bezier curve 
has self-intersection. For the other control polygon there are 
weights (e.g. Fig. |l(b)| i such that the resulting Bezier curve 



has a point of self-intersection. The cited works provide 
conditions which imply no self-intersection. Our purpose 
is different: We give conditions on the control points for 
2D patches which are equivalent to there being no self- 
intersection for any choice of positive weights. 

The basic units in the geometric modeling of sur- 
faces are rational Bezier simplices and tensor product 
patches. Krasauskas [8| introduced toric Bezier patches as 
a natural extension of classical rational patches and their 
higher-dimensional generalizations, the Bezier simploids by 
DeRose, et al. [4|. The theory of toric patches is based 
upon real toric varieties from algebraic geometry (9), and 
they provide a general framework in which to pose many 
questions concerning classical rational patches. 

To study dynamical systems arising from chemical reac- 
tion networks, Craciun et al. [3| prove an injectivity theorem 
for certain maps. This was adapted in [2] to give a geometric 
condition on a set of control points which implies that the 
resulting toric Bezier patch has no self-intersection, for any 
choice of positive weights. That result contains a minor flaw 
in that it only guarantees injectivity in the interior of a patch. 
We correct that flaw, at least for 2D patches, showing that the 
condition from [2| plus the mild additional hypothesis that 
the vertices correspond to distinct control points is equivalent 
to injectivity for every choice of positive weights. 

In Section 2, we introduce toric Bezier patches as gen- 
eralizations of the classical rational patches. In Section 3 
we explain our condition and sketch its equivalence to the 
injectivity of every 2D patch with a given set of control 
points, for all possible weights. More details, including 
examples of the geometric arguments of Lemma 13.51 and 
Corollaries 13 .6 1 and 13 . 7 1 will be added in the complete version 
of this paper. We conclude some remarks on how to check 
this condition, argue that it is in fact quite natural, and 
interpret it in terms of piecewise linear maps. 

While our main interest is in establishing a criteria valid 
in 3D, and in fact in all dimensions, we currently do not 
know how to add hypotheses to the condition of [2| so that 
the result will be equivalent to injectivity for any choice of 
weights in 3D. 






Figure 1. Cubic Bezier curves. 



II. Toric Bezier patches 

Let A C 1? be any finite set of integer lattice points. 
Its convex hull A A is a polygon whose vertices are lattice 
points. This polygon is also defined by its edge inequalities, 

&a = {{x, y) G R 2 | < hi(x, y),i = l,...,£}, 

where hi(x) = a{X + b t y + a are linear polynomials with 
integer coefficients and (a,i,bi) is relatively prime. 

For each integer lattice point a G A, Krasauskas J8) 
defined the toric Bernstein polynomial 

/3a(a?) := h 1 (x) h ^ ) h 2 (x) h2 ^ ■ ■ ■ h e (x) ht{& \ (1) 

These toric Bernstein polynomials are non-negative on A .4, 
and the collection of all /3 a has no common zeroes in A^. 

Let R^ be R^ 1 with coordinates (w a G R> | a G A) 
indexed by elements of A. 

Definition 2.1: Let A C 1? be a finite set. A toric Bezier 
patch associated with A requires an assignment / : A — > R d 
(d = 2, 3) of control points and a choice of weights w G R> . 
The toric Bezier patch F w : A_a — > R d is the function 



F w (x) = F AJtW (x) 



(2) 



written F w as A and / are understood. 

The degree of a toric Bezier patch is encoded in its 
domain, differing from the classical patches as developed 
in J5) . These two types of patches share many properties, 
which is explained in J8], J9). Two properties in particular 
are important for us. 

One is the convex hull property, that the image of A^ 
under F w is contained in the convex hull of the control points 
f(A) with F w (b) = /(b) if b is a vertex of A^, and the 
other is the boundary property, that the restriction of F w 
to an edge 5 of A_4 is a rational Bezier curve, defined by 
control points and weights corresponding to lattice points of 
5. 

The boundary property may be seen directly by consider- 
ing the restriction to an edge. For the convex hull property, 
note that as w & f3 & (x) is nonnegative, F w (x) is a convex 
combination of the control points, and if b is a vertex, 
then /3 a (b) is zero unless a = b. Since the toric Bernstein 
polynomials are strictly positive on the interior of A (and 



those corresponding to an edge S are strictly positive on the 
interior of 8), we may deduce a little more. 

Proposition 2.2: The image of the interior of A lies 
strictly in the interior of the convex hull of the control points 
f(A), and the image of the interior of an edge S lies strictly 
within the interior of the convex hull of f(S n A). 

Toric Bezier patches include the classical Bezier patches 
and some multi-sided patches such as Warren's polygonal 
surface ifTUI which is a reparameterized toric Bezier surface. 

Example 2.3 (Tensor product patches): Let m,n be pos- 
itive integers. Let A be the integer points in the m x n 
rectangle A := ■ < i < m, < j < n}. Then the 

corresponding toric Bernstein polynomials (Q~|) are 

P (l , 3) (x,y) := a: i (ro-zr~V(n-!O n ~ i J (3) 



and the toric Bezier patch (O (with weights 



J i,3 



% )0 



is the rational tensor product Bezier patch of bidegree (m, n) 
after the simple reparameterization s = x/m, t = y/n. 

Example 2.4 (Triangular Bezier patches): Let m be a 
positive integer and A be the integer points in the triangle 
with vertices (0,0), (m, 0), and (0, m), A := \ 
< < in — i, j}. The corresponding Bernstein 

polynomials ([T]) are 

Pi,j(x,y) = xty* (m - x - y)™-*-' . 

Then the toric Bezier patch (|2]) (with weights wtj = 
n .,, m! - — ^t) is the rational Bezier triangle of degree m after 

i\jl(m — i— j)\ ' fc & 

the simple reparameterization s = x/m, t = y/m. 

III. INJECTIVITY OF 2D TORIC BEZIER PATCHES 

Given a finite set A G 1? and a choice /: A — > R 2 of 
control points, we consider the injectivity of toric Bezier 
patches as mapping functions F w : A _4 i-> R 2 ©, for all 
choices w G R> of positive weights. 

Affinely independent points ao , ai , &2 determine an ori- 
entation via the ordered basis ai— ao, a.2 ao of R 2 . 

Definition 3.1: A choice /: A — > R 2 of control points is 
weakly compatible if 

1) There are affinely independent points ao,ai,a2 of A 
such that /(ao), / (ai), /(a2) is also affinely indepen- 
dent, and 

2) For any affinely independent points a , a.[ , sl 2 of 
A with the same orientation as ao,ai,a2, if 



/(ad), /(a'i), /( a 2 ) is a l so affinely independent, then 
it has the same orientation as f(&o), /(ai), f(&2)- 
Fig. 12 shows three sets of labeled points, indicating 
assignments between them. The assignment between the first 
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Figure 2. Weak compatibility. 

two sets is weakly compatible, but neither assignment to the 
third set is weakly compatible. 

We state Theorem 3.5 of [2| for R 2 , which is their main 
result on injectivity of toric Bezier functions (it holds in any 
dimension). Write for the interior of A_^. 

Theorem 3.2: The map F w : H> R 2 is injective for 
all w G R> if and only if the assignment /: A — > R 2 is 
weakly compatible. 

In Q, the authors incorrectly stated this result as F w is 
injective on all of A^, even though their proof was only 
valid for the interior of the convex hull. Their proof showed 
that F w has no critical points in the interior, which shows 
that it is an open map on A^. 

This is the best possible result with these hypotheses: 
Consider a bilinear patch where two control points coincide. 
Specifically, let A = {(0, 0), (0, 1), (1, 0), (1, 1)} and sup- 
pose that the control points are {(0, 0), (0, 1), (1, 0)}, where 
/(a) = a, except that /(l, 1) = (1,0). This assignment of 
control points is weakly compatible, but F w collapses the 
edge between (1,0) and (1,1) to the point (1,0). 
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This example shows that more hypotheses are needed to 
ensure that F w is injective on A_4, and those hypotheses 
should imply that faces of A_4 are not collapsed. In fact, 
this is the only additional hypothesis needed. 

Definition 3.3: A choice / : A — > R 2 of control points is 
compatible if it is weakly compatible, and no two vertices 
have the same image under /. 

We state our main result. 

Theorem 3.4: The map F w : A^ h-> R 2 is injective for 
all w G R> if and only if the assignment /: A —> R 2 is 
compatible. 

If a <E A is a vertex of A_4, then F w (sl) = /(a). 
Theorem 13.21 together with this observation, shows that 
if F w is injective for all w G R^, then /: A — >• R 2 is 
compatible. 



For the other implication, suppose that /: A —> R 2 is 
compatible. We show that the assumption that F w is not 
injective leads to a contradiction. 

We first make several observations about the relative 
positions of the points /(a) for a G A which are im- 
plied by compatibility. Composing with a reflection of 
R 2 if necessary, we may assume that if arj,ai,a2 and 
/(ao), / (ai), /(a 2 ) are both affinely independent, then they 
induce the same orientation on R 2 . 

Let 6 be an edge of A^. There is some triple of points 
d, d',a of A with /(d), f(d'), /(a) affinely independent 
where d, d' G S and a g" S. Indeed, if there are no such 
triples, then every point of f(A) lies on every line segment 
between two distinct points of f(6 HA), which implies that 
the points of f(S n A) are collinear and the line they span 
contains f(A), which contradicts the first condition for weak 
compatibility of Definition 13.11 This argument requires that 
there be at least two distinct points of f(8 n A), which 
follows as the endpoints of S (which are vertices of AX) 
are mapped to different points under /. 

Suppose that we list the points do, di, . . . , d m of 6 P\ A 
so that if a G A \ S, and i < j, then dj, dj, a are positively 
oriented. Then either / (dj), /(dj), / (a) are collinear or pos- 
itively oriented. Since there must be at least one such triple 
with /(dj), /(dj), /(a) affinely independent, we deduce the 
following. 

Lemma 3.5: Every control point f(A \ S) lies in the 
intersection of closed halfspaces 

{x G R 2 | f(di),f(dj),x are positively oriented} 

for i < j with /(dj) ^ /(dj), and this intersection has a 
nonempty relative interior. 

Corollary 3.6: For every edge S of A_4 and every b G 
A \ 6, the control point /(b) does not lie in the relative 
interior of the convex hull of f(5 n A). 

To see this, note that the intersection of halfspaces of 
Lemma 13.51 is either interior or exterior to the convex hull 
of f(5 n A), and if it is exterior, then it is separated from 
the relative interior of the convex hull by a line. If there is 
an edge S so that this intersection lies in the interior of the 
convex hull of f(5P\A), let 6' be a different edge. Then the 
positions of the points of S n A relative to the intersection 
of halfspaces for 5' leads to a contradiction. 

Corollary 3. 7: If / : A —5- R 2 is compatible, then the 
restriction of F w to any edge S of A^ is injective. 

To see this, fix an edge S and consider the intersection of 
halfspaces of Lemma [33] This intersection is exterior to the 
convex hull of f(6<l A) and so consists of an unbounded 
polyhedron, P. Consider the orthogonal projection ir : R 2 — >■ 
R along an unbounded direction of P. Then the map ir o 
/:<5n.A— ^Risa weakly compatible choice of control 
points for 6 f) A, and so the map it o F w restricted to the 
edge 5 is injective, by Theorem 13.21 But this implies that 
the restriction of F w to 6 is injective. 



Proof of Theorem \3.4\ We suppose that /: — 5- IR 2 is 
compatible and that F w is not injective. Let x, y € A ,4 be 
distinct points with F w (x) = F w (y). 

First, neither x nor y can be a point of A^. To see this, 
suppose that x G and let V be a neighborhood of x 
in A^4 whose closure does not contain y. Then F W (V) is 
an open set containing F w (x) = F w {y), so F^Q/) \ V 
contains an open subset U of y in A .4. But then points of 
U fl A^ are mapped by F w to points of F W (V), and so 
F w is not injective on the interior of A, which contradicts 
Theorem 13.21 as the choice / of control points is weakly 
compatible. 

Thus x and y are points of some edges of A ,4. They 
cannot be points of the same edge 5, for then the restriction 
of F w to 6 is not injective, contradicting Corollary 13.71 Thus 
they are points of different edges, x £ S and y £ 8' with 
5 ^ 5'. We cannot have one of them be an interior point of 
its edge, for then the relative interiors of the convex hulls of 
f(S ("1 A) and f(S' n A) meet, contradicting Corollary [3~6l 

The only possibility left is that x and y are vertices of 
A, but then F w (x) = f(x) and F w (y) = f(y), which are 
different, as the choice / was compatible. ■ 

Remark 3.8: By definition, to check weak compatibility 
for 2D patches, it suffices to check determinants for each 
triple of points of A and the corresponding control points, 
giving a simple (#(_4)) 3 algorithm. The complexity may be 
reduced if we start from a triangulation of A_4, or with care- 
ful bookkeeping. Such triangulations can be obtained from 
control nets for tensor product patches or Bezier triangles. 
We will treat the complexity of checking weak compatibility 
in the complete version of this extended abstract. 

Mapping functions that are weakly compatible exist; for 
example the identity assignment of control points is weakly 
compatible. A designer may choose weakly compatible 
control points for aesthetic or other reasons. For example, 
if only a few control points are moved such as in image 
warping, morphing, or reparameterization, then the control 
points may be weakly compatible by design, or else only a 
few determinants need to be computed. 

For any triangulation of A, the assignment of control 
points induces a piecewise linear map to the image. This 
piecewise linear map is injective (except possibly collapsing 
an interior simplex) for every such triangulation if and only 
if the assignment of control points is weakly compatible. 

IV. Conclusions 

In this paper, we study the injectivity of toric Bezier 
patch geometrically. We present a simple condition on a set 
of control points which implies that the resulting 2D toric 
Bezier patch is injective, for any choice of positive weights. 
For higher dimension, the best result remains Theorem 13.21 
by Craciun et al. in [2| (Theorem 3.5 in 121). We plan to 
continue this investigation of injectivity for 3D and higher 
dimensions in a future publication. 



Acknowledgements 

The authors thanks to Tim Goodman and Keith Unsworth 
for providing their paper |7 j. Research of Sottile is sup- 
ported in part by NSF grant DMS-1001615, and the Institut 
Mittag-Leffler, Djursholm, Sweden. Research of Zhu is 
supported by the NSFC (Grant Nos. 10801024, 11071031, 
and U0935004), the Fundamental Research Funds for the 
Central Universities (DUT10ZD112, DUT11LK34), and 
the National Engineering Research Center of Digital Life, 
Guangzhou 510006, China. 

References 

[I] Y. Choi and S. Lee, "Injectivity conditions of 2D and 3D 
uniform cubic B-Spline functions," Graphical Models, vol. 62, 
2000, pp. 411-427. 

[2] G. Craciun, L. Garcfa-Puente and F. Sottile, "Some geometrical 
aspects of control points for toric patches," in: Mathematical 
Methods for Curves and Surfaces, M. Dahlen et al. Eds, 
Lecture Notes in Computer Science, vol. 5862, Springer, 
Heidelberg, 2008, pp. 111-135. 

[3] G. Craciun and M. Feinberg, "Multiple equilibria in complex 
chemical reaction networks I. The injectivity property," SIAM 
Journal on Applied Mathemathics, vol. 65, 2005, pp. 1526- 
1546. 

[4] T. DeRose, R. Goldman, H. Hagen, and S. Mann, "Functional 
composition algorithms via blossoming," ACM Transactions 
on Graphics, vol. 12, 1993, pp. 113C135. 

[5] G. Farin, Curves and Surfaces for Computer Aided Geometric 
Design, Computer Science and Scientific Computing, San 
Diego: Academic Press Inc., 1997. 

[6] M. Floater, "One-to-one piecewise linear mappings over tri- 
angulations," Mathematics of Computation, vol. 72, 2003, pp. 
685-696. 

[7] T. Goodman and K. Unsworth, "Injective bivariate maps," 
Annals of Numerical Mathematics, vol. 3, 1996, pp. 91-104. 

[8] R. Krasauskas, "Toric surface patches," Advances in Compu- 
tational Mathematics, vol. 17, 2002, pp. 89-113. 

[9] F. Sottile, "Toric ideals, real toric varieties, and the moment 
map," in: Topics in Algebraic Geometry and Geometric Model- 
ing, Contemp. Math., vol. 334, Amer. Math. Soc, Providence, 
RI, 2003, pp. 225-240. 

[10] J. Warren, "Creating multisided rational Bezier surfaces using 
base points," ACM Transactions on Graphics, vol. 11, 1992, 
pp. 127-139. 

[II] G. Wolberg, Digital Image Warping, Los Alamitos: IEEE 
Computer Society Press, CA, 1990. 



